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Reduction of convective beat transfer with altitude dictates unusually large heat exchangers for piston- 
engined high-altitude aircraft. The relatively large aircraft drag ftnction associated with cooling at high 
altitudes makes the efficient design of the entire heat exchanger installation an essential part of the 
aircraft’s aerodynamic design. The parameters that directly influence cooling drag are developed to the 
context of high-altitude flight Candidate wing airfoils that incorporate heat exchangers are ««“““• 
Such integrated wing-airfoil/heat-exchanger installations appear to be attractive alternatives to isolated 
heat-exchanger installations. Examples are drawn from integrated installations on existing or planned 
high-altitude aircraft. 


Nomenclature 


A r > A h 

= radiator frontal area and heat 
transfer area 

C f 

= radiator skin friction coefficient 

Cl, Cp 

= profile lift and drag coefficients 

Cp 

= pressure coefficient 

c 

= wing chord 

Cp 

= specific heat at constant pressure 

£>co« 

= radiator core drag 

H 

= shape parameter, S*/6 


= heat rejection rate 

h 

= radiator height, enthalpy 

K p K h 

= scaled friction and heat transfer 
coefficients 

k 

= thermal conductivity 

l 

= radiator air passage length 

M 

= Mach number 

m 

= radiator air mass flow, piVjA, 

P 

= radiator pressure-drop parameter, 
Ap/<? t 

Pr 

= Prandtl number, c p ix/k 

P 

- pressure 

q 

= dynamic pressure, t pV 2 

r h 

= hydraulic radius 

St 

= Stanton number 

T 

= temperature 

tt, V 

= boundary- layer velocity components 

V 

= velocity 

y 

= boundary -layer coordinates 

y 

= ratio of specific heats 

6*, e 

= displacement, momentum thickness 


= viscosity 

f 

= passage length parameter 

p 

= density 

a 

= radiator blockage factor 


Subscripts 

e = boundary-layer edge 

0 = stagnation 

1, 2, 3 = front, back, downstream of radiator 

oo = free stream 
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I. Introduction 

S UBSONIC ultrahigh-altitude aircraft are currently being 
developed for use in a variety of missions, notably for in 
situ atmospheric science measurements. 1 3 Although previous 
high-altitude aircraft have not flown above 70,000 ft (e.g., 
Boeing’s Condor 4 ), cruise altitudes of up to 30 km (100,000 
ft) are being considered for the new science mission-oriented 
aircraft. Practical constraints of low development and unit cost, 
and the fact that turbofan and turboprop engines begin to suffer 
large power losses near 20 km (65,000 ft), dictate the use of 
turbocharged reciprocating engines for such aircraft. A recip- 
rocating engine used at high altitudes must in practice be liquid 
cooled and the size of the radiator must be considerably larger 
at high altitude than at sea level. Russ et al. 5 have outlined the 
system integration issues involved in the development of such 
radiators. Careful attention must be paid to the aerodynamic 
design of the radiator installation, particularly for long-range 
aircraft whose performance is sensitive to the associated par- 
asite drag. . 

The large radiators typically required for high-altitude air- 
craft can adversely impact the flow around neighboring aero- 
dynamic components. For example, in the case of a wing lead- 
ing-edge inlet, an ineffective inlet geometry can have very 
adverse effects on the wing airfoil, as demonstrated in previous 
experimental studies. 6,7 This makes an isolated nacelle radiator 
installation less risky and more attractive. On the other hand, 
effective integration of the radiator with the wing airfoil also 
offers possibilities for favorable interactions and greater com- 
pactness, with possibly lower overall drag and lower system 
weight than with an equivalent isolated installation. 

This article examines the aerodynamics of closely coupled 
airfoil/radiator installations and develops criteria for selecting 
important installation and operating parameters. A computa- 
tional two-dimensional method for predicting performance of 
the entire airfoil/radiator configuration will be used to illustrate 
the pros and cons of several alternative configurations being 
employed on operational and planned ultrahigh-altitude 
aircraft. 

H. Characteristics of Low-Drag Heat 
Exchanger Installations 

The ideas behind low drag radiator installations have been 
known since before World War II (see Ref. 8). Hoerner dis- 
cusses installations employed to date. The basic principle is to 
decelerate the cooling airflow, pass it through the radiator at 
low speed, and then accelerate it back to ambient pressure as 
shown in Fig. 1. In practice, the final streamtube area A* is 
very close to the geometric nozzle area. 
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Fig. 1 Cowled radiator and radiator core section. 


The drag D «« of the radiator core is defined as the cooling 
airflow’s rate of change of momentum: 

D «« = - V 3 ) = p,V,A,(V. - V 3 ) (1) 

Adding to this will be the profile drag of the cowl D and 
perhaps some interference drag due to unfavorable interac- 
tions. In a well-designed heat exchanger system, the D ^ and 
D COVfl drag components should be predominant, and will be the 
primary topic here. 

A. Radiator Core Geometry 

In most liquid/air radiator cores, such as those found in au- 
tomobiles, the liquid coolant flows perpendicular to the airflow 
in flattened tubes that are separated by corrugated metal fins 
that form cooling-air passages, as shown in Fig. 1 . The cross- 
sectional size of the air passages is commonly characterized 
by r h , defined in terms of the ffontal/heat transfer area ratio, 
and the streamwise flow length / (i.e., the thickness) of the 
radiator: 


/ <J ^ r _ passage area 
A h passage perimeter 


( 2 ) 


This r h definition reduces to half the passage width for a long 
rectangular section and to one-quarter of the passage diameter 
for a circular cross section. The presence of the liquid-carrying 
tubes, represented by o-, results in the flow in the passage inlet 
being increased from V l to VJcr. For most liquid/air radiators, 
<r is close to unity. 

For a compact radiator it is desirable to increase the internal 
radiator wetted area A h by decreasing r h and/or increasing l. 
However, once the thermal and momentum boundary layers 
begin to merge, as shown in Fig. 1, further increases in heat 
transfer are smaller compared to the increases in skin friction 
and radiator pressure drop. This effect can be examined via 
the average friction coefficient and Stanton number: 


2f _ oA r 2A p 
P.CV'./O-) 2 = 17 Pi(Vy<r) 2 


St = 


c P (T r - r,)piV,/o- 


(3) 

(4) 




□ 



Fig. 2 Average and local Reynolds number-scaled heat transfer 
and skin friction coefficients vs scaled passage length parameter 
for a two-dimensional channel. 


so that for the limiting case of zero blockage (cr = 1), the net 
pressure drop across the radiator would be A p. 

Figure 2 shows local and average Reynolds number-scaled 
friction and heat transfer coefficients 



(5) 

( 6 ) 


vs the scaled £ 


>2 m 2. Mi 
* r h Pr p { V,r h 


(7) 


for laminar entrance flow in a two-dimensional channel. The 
curves were computed using a finite difference thin shear layer 
code with a small Mach number and Pr = 0.72, although the 
dependence on the Prandtl number is very weak. The limiting 
case f — ► 0 corresponds to laminar fiat plate flow, for which 


Kf=K k = 


0.664 

0.332 


average 

local 


with K f =K h constituting the Reynolds analogy. 

Past f — 1.0, which is comparable to the entrance-flow 
length for a two-dimensional channel, the friction begins to 
increase more rapidly, while the heat rejection begins to level 
off and eventually decrease. This then sets a practical upper 
limit for the length of the radiator and often is the limit of 
measured C f and St radiator core data, such as that presented 
in Kays and London. 10 Core data for real radiator geometries 
differ quantitatively somewhat from the curves in Fig. 2, but 
still have the same general behavior. 

The choice of f — C(l), which defines an entrance-flow 
radiator, is likely to be effective in a high-altitude aircraft 
where it is necessary to minimize the inevitably large radiator 
frontal area, but without incurring an excessive pressure drop. 
In any situation, reducing the radiator drag can be achieved 
with a smaller £ while obtaining a more compact radiator with 
a smaller frontal area can be achieved with a larger £ 


The average shear stress f is related to the pressure drop cor- 
rected for the increasing momentum flux of the developing 
velocity profile in the passage 

Ap = A (p + J pV 2 


B. Pressure Loss — Drag Relation 
For the case of low freestream Mach number and frictionless 
radiator, Rauscher 11 has explicitly related V 3 and the drag D ^ 
to the heat rejection and exit area. For fully compressible flow, 
the dependence of the radiator drag on all of the operating 
parameters cannot be explicitly written down, but is implicitly 
determined by the governing equations for compressible quasi- 
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one-dimensional channel flow with friction and heat addition, 
as discussed in McGeer et al. 12 and Russ et al. 5 These equations 
can be readily solved by Newton iteration for the flow varia- 
bles (p, p, V) ixy The numerical two-dimensional simulations 
presented later duplicate such fully compressible results. Here, 
Rauscher’ s low-speed treatment will be extended to nonzero 
friction, to clearly identify the relevant parameters that influ- 
ence the radiator drag. 

Rauscher assumes that the only significant density change 
occurs across the radiator, 


Pi “ Pi ~ P 3 


and the drag power and equivalent radiator drag coefficient 
can then be estimated as 



where the approximation pjp x =* p«/p« has been made. Equa- 
tions (14) and (15) clearly identify the pressure coefficient 
drop across the radiator 


and employs the incompressible Bernoulli equation upstream 
and downstream of the radiator. His momentum equation 
across the radiator is extended here to include the frictional 
pressure drop parameter P> which is experimentally measured 
for a given radiator core and Reynolds number. For numerical 
simulation purposes, P is separated into the core inlet P h core 
exit P„ friction P fy and acceleration terms, as suggested by the 
pressure-drop model of Kays and London 10 : 


( 8 ) 


Pi - Pi - %P\V\P 


The friction term is often dominant and is simply related to 
the average skin friction coefficient and blockage factor: 

P f =(C f (* 2 )(l/r H ) (10) 

Typical overall pressure drop parameter values for liquid/air 
radiator cores are in the range 5 < P < 20. 

In his frictionless model, Rauscher employed only the ac- 
celeration term 


P = 2[(V 2 /Vj) - 1] 

and gave explicit (and elaborate) expressions for V u V 3 , and 
Deo,, for a given radiator, exit area, flight conditions, and heat 
rejection. For design purposes, it is more useful to consider 
the radiator velocity ratio V x /V„ as a free parameter, since this 
controls the radiator air mass flow, and immediately gives the 
radiator velocity and density ratios: 


AC,, * C„ - C Pl - P(V x tV„f (16) 

as the decisive parameter in controlling radiator drag, since the 
mass flow m ~ A r V x is strongly constrained by the requirement 
to reject the required beat load. The second term, proportional 
to HMi, is the ramjet thrust, which may significantly offset 
the first friction drag term in some installations. 

C. Heat Transfer — Drag Relation 
Sizing of the radiator requires a quantitative relation be- 
tween the known heat load H and the other radiator parame- 
ters, and the mass flow m = PiV t A r in particular. The experi- 
mental data in Kays and London 10 indicate that for typical 
plate/fin radiators in the laminar flow regime, the skin friction 
and heat transfer coefficients closely follow the Reynolds num- 
ber-dependent scaling implied in Fig. 2. Combining Eqs. (2) 
and (4) gives a direct relation between heat rejection and mass 
flow: 


H = £K*c p (T r - T x )Pr~ M piV x A r (17) 

From Fig. 2, the scaled Stanton number takes on the values 
K h ~ 0.7 -0.9 for a two-dimensional passage, but in practice 
would be taken from measured data for a particular radiator 
core. 

The fully developed flow case £» 1, considered by Capon, 8 
has a clear loss in effectiveness (smaller K A ) due to the passage 
centerline temperature asymptoting to the radiator metal tem- 
perature 7V, at which point further heat transfer cannot occur. 
The T r itself can vary in the range between the incoming air 
temperature and the liquid coolant temperature 

7V < T r < T ^ (18) 


Yi = 2±~Ti^ ! + — l£i (id 

V, p 2 7V my p t 

The net pressure drop parameter P can then be estimated from 
the Kays and London model [Eq. (9)], or preferably taken from 
measured data if available. Combining the radiator momentum 
Eq. (8) with the inlet and outlet Bernoulli and continuity equa- 
tions finally gives the exit velocity 

(vyv .) 2 = (p,/p 2 )[ i - Pivjv.f] (12) 


depending on the coolant circulation rate. Of course, T, can 
also vary across the radiator if the liquid coolant is not cir- 
culated and mixed rapidly. These effects cannot be considered 
here, but can in general be lumped into a net effective K„. 

The entrance-flow radiator heat rejection relation [Eq. (17)], 
combined with the radiator drag definition [Eq. (1)], results in 
a direct relation between the heat rejection and drag power 


D^V. = 


HVlPr^ 
£K*c,(7V - 7j) 



(19) 


and the radiator drag then follows from its definition [Eq. (1)]. 

A useful further approximation is to consider the case of a 
small exit velocity defect and small density change. 


Using the small-defect approximation [Eq. (12)], this becomes 


[ viPr M P (v A 2 

De " V " “ H ~ T.) 2 (vj 



( 20 ) 


V 3 /V. =* 1, p,/p 2 =* 1 


in which case we have 


If the friction term in Eq. (9) is the dominant contribution to 
P, as is usually the case, the drag power can be further sim- 
plified to 


Vs 
V . 




, 7 - 1 H P- 
2 m yp x 


(13) 


D^V. = H 


VlPr m 1 A>/V,\ 2 
c,(7V - T.) <x 2 K, \vj 



( 21 ) 
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Fig. 3 U.S. Standard Atmosphere density and temperature profiles, corresponding radiator velocity ratio required for cooling regulation, 
and associated drag power vs altitude for three radiator temperatures and three minimum velocity ratios. Assumes ip. Vi = 480 Pa = 
10 Ih/ft 2 and K f /<r 2 K k = 10. 


Equation (21) quantifies the benefits of using a small velocity 
ratio V,/V«, of selecting a radiator core with a small K f /K h 
ratio and near-unity blockage factor o\ and of using a large 
radiator temperature T r . 

D. Cooling Regulation Considerations 
The discussion up to now has been aimed at the sizing of 
an ideal radiator for a particular altitude. In practice, the ra- 
diator is of one fixed size, but it must still reject heat at the 
same rate at all altitudes, such as in a typical full-power climb- 
out. Recasting Eq. (17) as 

H = KJ>r~™ c p (T, - r.) (^j p.VA (22) 

shows that for a given radiator geometry and heat load, it is 
necessary to hold 

(T, - r.Xp.VV,) 1 ' 2 ~ const (23) 

with altitude and dynamic pressure. 

Constraint (23) can be satisfied either by regulating T r (by 
varying the liquid coolant circulation rate), or by regulating 
the radiator velocity ratio Vj/V* (by varying the nozzle area 
A 3 with a cowl flap). Regulating the coolant flow is undesirable 
in practice, since at low circulation rates the engine may en- 
counter large thermal stresses from the returning excessively 
cold coolant, or the coolant may freeze and block the radiator. 
In contrast, regulating the velocity ratio permits maintaining 
the same ideal coolant temperature with varying altitude and 
is much more attractive for reliable engine operation. 

It is reasonable to assume that the conditions at station 1 are 
nearly at stagnation: 

PI - Pb = P«{l + [( 7 - l)/2 

Hence, constraint (23) implies that for proper cooling regula- 
tion it is necessary to vary the radiator velocity ratio as 

V, {! + [(?- i)/y](<7~/p~)r Iy/< " 1>} 

V. ~ (T r - r„u + Kr - \)iy](qJp~)})\M~) m ^ { } 

using yM 2 = 2 qlp and assuming that p, - T. A more accurate 
p.(7) relation could be used. Figure 3 shows the required ve- 


locity ratio defined by Eq. (24) relative to its minimum value 
as a function of altitude, for three radiator temperatures and at 
fixed dynamic pressure. The radiator temperature does not sig- 
nificantly impact the range of Vi/V« required for cooling reg- 
ulation, but it affects the heat rejection per frontal area. In all 
cases, flight near the tropopause has the largest available cool- 
ing and will likely be the design point for the smallest V,/V« 
tolerable by the heat exchanger installation. 

To determine the radiator drag variation with altitude, it is 
necessary to require that the velocity ratio V,/V* in Eq. (20) 
varies with altitude according to the regulation constraint (24). 
The approximate drag-power Eq. (21) can be expressed as 

^ (2 gjpjpr™ i_ k, fvX 

H c p (T r - T-{1 + l(y - l)/2 ]Ml}) a 2 K h [vj 


with Vt/V« determined by Eq. (24) once the minimum value 
( Vi /V»)mir> is chosen. Since Ti is typically comparable to me- 
chanical engine power, Eq. (25) gives the fractional power loss 
(or gain) due to radiator drag or thrust. This is plotted in Fig. 
3 for the three radiator temperatures and three minimum ve- 
locity ratios. The negative ramjet term in Eq. (25) is also 
shown. The curves were generated assuming 


— 10 
2 o- 2 K„ 

which is typical for liquid/air heat exchangers. Rather small 
minimum velocity ratios at the tropopause are required if the 
radiator drag power is to be kept at a reasonably small level 
at much higher altitudes. Note also that for such designs, a net 
thrust is obtained from the ramjet term -[(y - l)/2]Afi over 
much of the operating altitudes. The major drawback to small 
minimum velocity ratios YVV« is that A r scales more or less 
inversely with V, for a given heat rejection rate, as can be seen 
from Eq. (17). This implies proportionately larger cooling sys- 
tem weight and larger duct cowls that also extract a propor- 
tionately larger profile drag penalty. Clearly, a careful tradeoff 
between these design factors is required. 
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HI. Numerical Simulation 

The numerical simulations presented here were generated 
with an extension of the MSES viscous/inviscid two-dimen- 
sional multielement airfoil code. 13 It employs a finite volume 
Euler formulation discretized on an intrinsic streamline grid. 
The viscous layers are represented by a two-equation lagged- 
dissipation integral formulation. The two formulations are 
strongly coupled through the displacement thickness and 
solved simultaneously by a global Newton method. 

A- Inviscid Flow Model 

The radiator in the simulation is placed along a row of cells 
that span the gap between two airfoil elements. Equation (8) 
together with the P definition [Eq. (9)] replaces the streamwise 
momentum equation at each radiator cell. The normal-momen- 
tum equation is replaced by the requirement that the flow di- 
rection at the exit (station 2) is along a specified direction, 
typically normal to the radiator face. This models the flow- 
straightening effect of the radiator core passages. 

The heat addition is simulated by specifying the stagnation 
enthalpy all along the streamtubes behind the radiator to be 
increased by an amount A/to. 

- h 0i * Aho = Him (26) 

The heat rejection rate H is prescribed, while the radiator mass 
flow m is, in general, unknown, and is generated as part of the 
solution. The temperature jump 

Ar= {A*o + [(V? - VD/21K1/C,) (27) 

is therefore a result of the solution as well. If it is found that 
the downstream temperature exceeds the coolant temperature, 

i.e.. 


T 2 = T x + AT > 

then an unrealistically large H is being specified. 


B. Viscous Flow Model 

A very important effect that must be modeled in a numerical 
simulation is the radiator’s thinning effect on the surface 
boundary layers that pass through it. Mehta 14 describes this 
effect for screens used to suppress separation in diffusers. Fig- 
ure 4 shows the pressure field set up by the radiator and the 
incoming boundary layer. Near the wall where the velocity is 
small, little or no pressure drop can be supported across the 
radiator, resulting in the low downstream pressure p 2 being felt 
in front of the radiator and accelerating the upstream boundary 
layer. The usual thin shear-layer approximation dpfdy =* 0 is 
obviously strongly violated, since at the upstream radiator face 
(Pe “ PwwuVp** 2 * “ P/2, which is typically much greater than 
unity. However, the effect can be estimated by adding an as- 
sumed radiator-induced pressure gradient term to the usual thin 
shear-layer momentum equation: 


du 


du 


pu Tx + pV 7y = 


+ 9r + 
dx dy 


- Ap f «*->* 
pm\ 1 dx 

(28) 


The radiator-induced pressure field sketched in Fig. 4 is as- 
sumed to propagate upstream with the length scale S compa- 
rable to the boundary-layer thickness and to scale with the 
local dynamic pressure pu 2 . Since the resulting acceleration of 
the boundary layer is an inviscid mechanism, the standard ap- 
proximation of neglecting streamwise diffusion is retained. 

Equation (28) produces the following modified integral mo- 
mentum and kinetic energy shape parameter equations: 


— = f e ~ H £ €n “* " {H + l)P 7 e ' 


dx 


d_ 

dx 


j*-x,ys 


(29) 



Fig, 4 Velocity profiles and pressure field associated with a 
boundary layer passing through radiator. 


d 2cj 


C, / 2 H**\ d - 

- — + [H - \ ~ - 77 — ) — €n u, 
20 \ H* ) dx 

( 2 H**\ P 

[ H - X ~ 1 F )2 


JX -Xrys 


(30) 


Here, c,and c 4 are the wall friction and dissipation coefficients, 
and //, H*, and H** are the conventional, kinetic energy, and 
density thickness shape parameters, respectively, as discussed 
by Drela and Giles. 13 The effect of the new additional terms 
containing P is dramatic. Even assuming the small value P * 
1.5, a roughly threefold reduction in the momentum and dis- 
placement thicknesses across the radiator is predicted for tur- 
bulent flow, consistent with the experimental observations 
mentioned by Mehta. 14 This strong effect can be examined by 
neglecting the skin friction and any u e variations just upstream 
of the radiator, in which case Eq. (29) becomes 


d(€n *) « -(// + (3D 

Integrating this for the typical values H =* 1.5, P =* 10 shows 
that enormous decreases in the momentum thickness are 
possible: 



exp[-(tf + \)P] - 10" U 


(32) 


In practice, the momentum thickness is virtually eliminated at 
the radiator, regardless of the particular choice of the length 
scale 6 (=6* + 56 is assumed). In fact, the numerical imple- 
mentation requires that P be artificially limited to moderate 
values, e.g., P < 3, to avoid numerical difficulties with a nearly 
vanishing boundary-layer thickness. 

As a design consideration, this boundary-layer thinning be- 
havior has a very beneficial effect in that the radiator has a 
strong tendency to suck out any separated fluid ahead of it, 
equalizing the flow rate all across the radiator face. This makes 
the overall configuration behavior relatively insensitive to the 
detailed design of the interior duct contours. 


C. Numerical Model Validation 
Figure 5 shows the computed streamline grid and C p con- 
tours for the integrated airfoil/heat-exchanger geometry used 
on the Perseus Unmanned Air Vehicle. The computation cor- 
responds to one of a set of wind-tunnel tests aimed at mea- 
suring the performance of the overall configuration. 16 A wire 
screen with a measured P = 6.7 was employed to simulate the 
radiator on the tunnel model. The heat rejection of the actual 
Perseus radiator was not simulated, but this has little impact 
on the aerodynamics in this case. Figure 6 shows the computed 
surface C p distributions and compares the computed and mea- 
sured drag forces. The computation separates the overall drag 
(momentum defect) into contributions from the viscous wake 
and from the inviscid radiator core flow wake: 


C D — 4* Co 


(33) 
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Fig. 5 Computed streamline grid and C, contours for Perseus airfoil/heat-exchanger test model configuration. 


2 in 


Cxw = (V. - V) dm (35) 

faVooC J eontiow 


In the experiment, the total C D was measured by a wake rake 
spanning the entire viscous and core flow profile, and C Dam 
was deduced from the static and total pressure coefficients 
and measured immediately behind the radiator screen. 
With the assumptions that C p 0| = 1 in front of the radiator, and 
that no further total pressure loss occurs in the core flow down- 
stream of the radiator, Eq. (1) reduces to 


C zw = 2 QUc)y 


- CM - Vcl) 


where hie is the radiator height to reference chord ratio. Figure 
6 shows that both the total C D and CtW are accurately pre- 
dicted. The rather low Reynolds number Re - 3 X 10 5 causes 
substantial transitional separation bubbles to appear on all sur- 
faces. The single static pressure measurement behind the ra- 
diator screen matches the computation very well, indicating a 
correct prediction of the radiator air mass flow. Figure 7 com- 
pares computed polars for Re = 3 X 10 5 and 2.5 X 10 5 with 
measurements at three angles of attack. The agreement is very 
satisfactory, giving confidence to using the present numerical 
simulation method for investigating design issues in integrated 
airfoil/heat-exchanger configurations. 


IV. Integrated Airfoil/Heat-Exchanger Installations 

Although numerous types of radiator installations have been 
employed on piston-driven aircraft, on high-altitude aircraft the 


options have been more limited due to the sheer size of the 
radiators. The Boeing Condor 4 and the Strato-2C 3 have their 
radiators contained entirely in large nacelles mounted on the 
wing. The Perseus UAV has its radiators inside a cowl system 
integrated with the aft portion of the wing airfoil, as shown 
earlier. This design was a retrofit into an existing wing struc- 
ture and the minimal impact of the aft-mounted installation on 
the wing airfoil was a strong constraint. The Theseus UAV 
currently in development will employ an alternative integrated 
installation, with the cowl and radiator near the leading edge. 
Here, the presence of the radiator and cowl requires significant 
redesign of the local airfoil. The Perseus and Theseus aircraft 
are designed for ceilings in excess of 25 km, and the isolated 
nacelles required to house their radiators would have been im- 
practically large and incompatible with the payload and engine 
placement. 

The aerodynamic advantage of an integrated airfoil/heat -ex- 
changer design lies in the fact that the apparent freestream 
velocity seen by the cowl/radiator system can be reduced by 
the presence of the wing airfoil. If the cowl/radiator configu- 
ration shown in Fig. 1 is placed on the bottom of the wing 
airfoil, then the V*. in Fig. 1 is reduced to something less than 
the true freestream velocity by the airfoil circulation. The cowl 
then experiences lower dynamic pressures, and hence, lower 
profile drag. Also, the amount of flow deceleration necessary 
to reach the required radiator velocity V x is smaller, allowing 
the cowl chord to be reduced, leading to further weight and 
drag reductions. On the other hand, integrated airfoil/heat-ex- 
changer installations can give rise to unintended boundary- 
layer separation that must be controlled with careful design. 

A. Aft- Mounted Underwing Cowl 

Locating the heat exchanger at the airfoil aft bottom surface 
as with the Perseus configuration shown in Fig. 6 is structur- 
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Fig. 6 Computed C, distributions and comparison of computed 
and experimental C D for Perseus airfoil/heat-exchanger test model 
configuration. 



Fig. 7 Comparison of computed drag polars with measurements 
for test model of Perseus wing-airfoil/heat-exchanger configura- 
tion. 


ally unintnisive and is attractive from potential-flow consid- 
erations. The local flow direction is essentially independent of 
the aircraft’s angle of attack, and the local effective freestream 
velocity seen by the cowl can be as small as 0.8V« for aft- 
loaded airfoils. The biggest source of difficulty is the wing 
bottom surface boundary layer, which cannot possibly nego- 
tiate the deceleration to the near-stagnation pressure at the ra- 
diator face. The suckdown effect examined earlier alleviates 
this to a significant extent, but its range extends only over a 
few boundary-layer thicknesses, and hence cannot suppress 
separation well upstream of the radiator. An effective solution 
is a splitter surface, which allows the wing boundary layer to 
bypass most of the cowl’s pressure field. 

An effective solution is a splitter surface, which allows the 
wing boundary layer to bypass most of the cowl’s pressure 
field. Numerical investigations indicate that the length and po- 
sitioning of the splitter is important. Figure 8 shows the inte- 
gral heat exchanger installation of the Perseus aircraft with 
three splitter length versions, together with computed displace- 
ment bodies and several velocity profiles. The short splitter 
doesn’t sufficiently protect the wing boundary layer from sep- 
aration, whereas the long splitter has a thicker boundary layer, 
which itself separates. The intermediate splitter length strikes 
a reasonable balance between separation on the wing and split- 
ter, and on this basis was selected for Perseus. 

At the low Reynolds numbers being considered, the largest 
drawback of the aft-mounted heat exchanger installation ap- 


pears to be the difficulty of obtaining a radiator velocity ratio 
V,/V« much below 0.2 without incurring significant separation 
either on the airfoil or the splitter. As a result, the overall drag 
level of the airfoil/heat-exchanger combination shown in Figs. 
6 and 7 is three to four times that expected of the clean airfoil 
at the same operating conditions. 

B. Front-Mounted Underwing Cowl 

Most of the problems inherent in the aft-mounted radiator 
can be largely overcome with the front-mounted installation, 
which does not subject any thick boundary layer to the cowl/ 
radiator’s pressure field and allows much smaller velocity ra- 
tios to be used. Figure 9 shows the airfoil/cowl/radiator con- 
figuration being developed for the large unmanned Theseus 
aircraft, which has a 27-km design ceiling and is being targeted 
for atmospheric science missions. This configuration is a log- 
ical extension of the leading-edge radiator intake commonly 





Fig. 8 Displacement bodies (shaded) and velocity profiles for 
Perseus airfoil/heat-exchanger installations with various splitter 
surface lengths, at typical near-ceiling flight conditions: a = 5 deg, 
Re e = 5 x 10 5 , Af. = 0.4, tilpV'A, = 0.5. thlpVA, = a) 0.242, b) 
0.240, and c) 0.223. 
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Fig. 9 Computed C, distributions for Theseus airfoil/beat-ex* 
changer configuration. 
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Fig. 10 Captured stream tube and surface C, variation for V,/V. 
— 0.05, 0.10, and 0.20 adjusted via cowl flap position. 


used on WW II aircraft, but with roughly four times larger 
radiator frontal area. Since most of the flow deceleration oc- 
curs ahead of the intake, there is little problem with flow sep- 
aration, and velocity ratios V x /V„ less than 0.1 can be readily 
obtained. In the calculation in Fig. 9, which corresponds to 
operation near the tropopause, the radiator core drag is virtu- 
ally nil even at the rather high pressure drop coefficient P = 
10, most of the frictional drag being cancelled by the ramjet 
effect. 

With the large variations in the velocity ratio V x /VL required 
between the tropopause and ceiling (Fig. 3), the detailed design 
of the duct intake lip shapes is quite important. Figure 10 
shows the behavior of the captured streamtube and the corre- 
sponding migration of the stagnation points as V X JV » =* ml 
p*>V*A r is varied from 0.05 to 0.20 by adjusting the exit area 
via the cowl flap. Most of the surface pressure variation occurs 
on the cowl, which is favorable in that of the overall 
section is only marginally affected. The presence of the cowl 
does reduce the relative to the clean airfoil, however, 
primarily due to the cowl’s negative camber. A possible com- 
plication in the design of the cowl is the rather low cowl- 
chord Reynolds numbers that can occur at ceiling. Fortunately, 
here the velocity ratio is typically set at its maximum and the 
pressure gradients on the cowl surface are most benign. A 
relatively open cowl flap will also incur the smallest pen- 
alty on the airfoil, which is also favorable with regard to ceil- 
ing performance. 

Wind-tunnel tests were performed on the Theseus front- 
mounted airfoil/radiator configuration to gain confidence that 
low velocity ratios could be attained at small chord Reynolds 
numbers. The cowl flap on the 10-in. chord model was hinged, 
allowing adjustment of the velocity ratios as would occur on 
an actual installation. A wire screen was again used to simulate 
the radiator. Figure 1 1 shows three drag polars corresponding 
to 0-, 5-, and 10-deg flap settings, producing V x /V . = 0.11, 
0.14, and 0.18 near the design condition of C L =* 1.1. Each 
polar is at the same chord Reynolds number of Re c = 4.5 X 
10 5 , roughly corresponding to a 27 km altitude. Also shown 
are the computed and measured (via stethoscope) upper-sur- 
face transition locations. These tests were of a preliminary na- 
ture, using a hand-built airfoil. The model airfoil differed 
somewhat from the design geometry, but the computations 
used the measured model coordinates for a direct comparison. 
A simple wire/strain-gauge balance used for the lift measure- 
ment and a total-pressure wake rake were used to measure the 
overall momentum defect. The experimental and computed C D 
curves in Fig. 1 1 actually represent the momentum thickness 
at x/c - 1.5, where the wake rake was located, to permit a 
direct comparison. The far-downstream momentum defect, and 
hence the true C Dl is computed to be about 10% smaller. 
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Fig. 11 Computed and experimental drag polars and upper-sur- 
face transition location of Theseus front-mounted radiator test 
configuration for 0-, 5-, and 10-deg cowl flap settings. 


These modest experimental results clearly indicate that 
rather small velocity ratios are viable on low Reynolds number 
installations. The predicted sudden C L break at stall is not con- 
sistent with the data, but the rounding off of the drag polar 
near is typical of hand-built models, which are prone to 
premature leading-edge C p spikes and consequent premature 
transition and separation at high angle of attack. 

At the smallest velocity ratio VJV<* = 0.11, at 0-deg flap, 
some separation is incurred on the cowl, especially at the lower 
a values where a significant suction spike occurs on the cowl 
leading edge. Figure 12 shows the cowl separation, which 
matches the separation observed in the experiment using a 
stethoscope. The separation adds significantly to the viscous 
profile drag component and this mostly cancels the ad- 
vantage of the smaller radiator drag component Ccw How- 
ever, this situation is not relevant for most applications, since 
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Fig. 12 Cowl separation induced by small velocity ratio V x fV m 
(dosed cowl flap) and small Reynolds number. 


such a small V,/V. would only be used at lower altitudes and 
much higher Reynolds numbers where the cowl is much more 
resistant to separation, such as in the tropopause-altitude case 
shown in Fig. 9. These test results indicate the limit of using 
small velocity ratios at small Reynolds numbers. 

V. Conclusions 

This article has discussed the issues in the design of low- 
drag heat exchanger installations for high-altitude aircraft. Pa- 
rameters that characterize the effectiveness of such installa- 
tions have been derived. The possibility of large drag penalties 
from poor installation designs has been shown to increase rap- 
idly with altitude above the tropopause. The requirements of 
heat rejection regulation via variable geometry have also been 
examined. 

A viscous/inviscid numerical simulation method for inte- 
grated airfoil/radiator flows has been demonstrated to be quite 
accurate at low chord Reynolds numbers and has proven useful 
in the development of candidate designs. The pros and cons 
of aft-mounted and front-mounted airfoil/heat-exchanger con- 
figurations have been examined with both computational re- 
sults and some preliminary experimental data. The front- 
mounted configurations appear to be inherently superior. All 
integrated installations demand close attention to proper de- 
sign. 
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